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Abstract: The current article completes our investigation of the hard-particle interac¬ 
tion by determining their distribution functionals. Beginning with a short review of the 
perturbation expansion of the free-energy functional, we derive two representations of 
the correlation functionals in rooted and unrooted Mayer diagrams, which are related by 
a functional derivative. This map allows to transfer the mathematical methods, devel¬ 
oped previously for unrooted diagrams, to the current representation in rooted graphs. 
Translating then the Mayer to Ree-Hoover diagrams and determining their automor¬ 
phism groups, yields the generic functional for all r-particle distributions. From this 
we derive the examples of 2- and 3-particle correlations up to four intersection centers 
and show that already the leading order reproduces the Wertheim, Thiele, Baxter solu¬ 
tion for the contact probability of spheres. Another calculation shows the failure of the 
Kirkwood superposition approximation for any r-particle correlation. 

Key words, integral geometry, fundamental measure theory, correlation functionals 


1. Introduction 

Density functional theory (DFT) for classical particles introduces direct correlation and 
distribution functionals. These two classes form the building blocks not only for the 
grand canonical potential and its perturbative expansion but also for the Enskog in¬ 
tegrals in non-equilibrium thermodynamics and the background potential in quantum 
mechanical calculations ITTIISI ITtI . But contrary to the correlation functions, which can 
be derived from molecular dynamic or Monte Carlo simulations, the correlation func¬ 
tionals are not readily available by numerical methods, which explains why so much 
less is known about this essential part of classical DFT. 

The two types of correlations derive from the free energy and the grand canon¬ 
ical potential respectively and are therefore canonically conjugate variables of their 
Legendre-dual potentials. This connection allows, at least in principal, to write any 
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thermodynamic expression in only one class of correlations and to use integral equa¬ 
tions to substitute the other ini. But solving these identities is a difficult task as the 
Ornstein-Zernike equation shows, whose only known analytic solution has been ob¬ 
tained by Wertheim, Thiele, and Baxter for the interaction of hard spheres ifTl ISSllTSl 
[^ . But despite its accuracy for low to medium packing fractions and its extension to 
the mean spherical approximation Cll, its solution is of limited use for more general 
geometries as the underlying mathematical methods do not generalize to non-spherical 
interactions. Another and even more fundamental restriction for the use of integral equa¬ 
tions is their dependence on the correlation functions instead of the functionals, which 
excludes their application for inhomogeneous fluids. 

To gain some insight into the structure of the correlation functionals, the current ar¬ 
ticle follows a different strategy and starts from the virial expansion in rooted Mayer 
diagrams mil. For a general potential this ansatz with its infinite number of integrals 
is inaccessible. But for the hard-particle interaction it has recently been shown that the 
free-energy functional expands in a fast converging series of intersection kernels, gen¬ 
eralizing previous results from Rosenfeld’s fundamental measure theory 11^124112511261 
ED, the transformation from Mayer to Ree-Hoover diagrams Ill81ll91l20ll . Wertheim’s 
derivation of the third virial coefficient II37II381I39II40I . and the Blaschke, Santalo, Chern 
equation from integral geometry Ill4ll5ll^ l^ . Its leading contribution coincides with 
the Rosenfeld functional, which has been shown to accurately predict the phase diagram 
of spheres and polyhedrons from low to medium packing fractions ll3ll9l fT0lfT4ll 15 II 1 61 
l27ll28lE0l . It is therefore a natural step to generalize these methods to the more general 
correlations functionals. 

As the direct correlation functionals derive from the free energy, they will be ignored 
in the following discussion. The same applies to the intersection kernels, whose deriva¬ 
tion is independent of the virial expansion. This leaves us to transfer the diagrammatic 
methods from the unrooted Ree-Hoover graphs to its rooted form, which will be done 
in Sec. where the generic distribution functional is derived. Examples will then be 
presented in Sec. where the 2- and 3-particle distributions for up to four intersection 
centers are given. It is then shown that already the leading term of the contact proba¬ 
bility between equally sized spheres agrees with the solution from Wertheim, Thiele, 
and Baxter. We conclude with a final comment on the applicability of the Kirkwood 
superposition approximation. 


2. The Generic Distribution Functional for Hard Particles 

The main difference between the free-energy and the distribution functionals is their re¬ 
spective virial expansion in Mayer clusters. Once this is known, it is a mere technicality 
to identify their intersection classes and to write the functional in intersection kernels. 
In the following, we will therefore first derive the Mayer representation of the r-particle 
distribution functionals p,j,or more conveniently of their normalized form g,,. 
translate them into Ree-Hoover diagrams, from which follows the generic correlation 
functional. 

A convenient starting point for the derivation of the Mayer representation is the per¬ 
turbation expansion of the grand canonical potential o. Introducing the hard and 
soft 0® contributions of the interaction potential — (j)^ + results in a correspond¬ 
ing splitting of the Boltzmann functions 

eij = efj+efjf^j = efj+XFij for A = l, 


(1) 
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where we introduced the auxiliary variable X for counting the number of F-terms. Next 
observe that the partition function of N particles is a fully e-bonded cluster integral of 
N labeled nodes iAr(e) and that the expansion of the product + ^Fij) yields a sum 
of products, with a subset of e-bonds replaced by F-bonds. Using the invariance of the 
partition function under relabeling of particle numbers, it is always possible to define 
a unique, labeled subgraph of r nodes and counting index k, such that each 

node is linked to at least one F-bond. The Taylor expansion of the partition function in 
X can now be written as a functional derivative of cluster diagrams 


~ X^FM ar,k ^ , n 5 

17^, II ’F) 


r=2k=\ 


^r,k\ 


SFrie) 


A=0 


( 2 ) 


where denotes the number of F-bonds, |f7| the total number of diagrams, and 
Gnk the number of inequivalent particle labelings. 

The Taylor expansion of the partition function up to second order includes the 
three leading diagrams, whose functional derivatives 


Zn(X). 


N=0^^ ■ i,;=l A'=0^''- 

= Eo[l+xJ + Y J + 


(3) 


‘ 2 ‘ 3‘4 


■ff{X^)] 


introduces the grand-canonical r-particle distribution functionals pP . Expanding its 
logarithm and setting A = 1 reproduces the well known perturbation expansion of the 
grand-canonical potential ini 


/3f2 = -^J dJhh -\J P^lhhfhhfhh ^yhhh 

_ 1 / (nF . . -nF nF f.S. dv ■ ■ ■ - 

g / ri\i2r 12,14)J i]_i2P 


(4) 


where an implicit sum over paired indices is understood. Higher order corrections are 
determined likewise by successive insertion of further F-bonds into the cluster dia¬ 
grams. 

The operator Dx provides a compact notation for the correlation functionals, where 
each diagram corresponds to exactly one pP ; . To change the representation 

into Mayer diagrams, one only has to insert F = e^f^, F I and to expand the 

graph in the virial series 




(5) 


n>r 


which yields a corresponding representation in terms of r-rooted Mayer diagrams FJ 


('■) 
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Fig. 1. The fourth virial order of the 2-particle correlation function in a) Mayer and b) Ree-Hoover diagrams. 
The continuous lines denote f-bonds, while the dashed lines correspond to e-bonds. The rooted points are 
always mutually e-bonded. 


where we omitted the hard-particle index and introduced the symmetry factor 
transforming from labeled to unlabeled graphs 022II34II . 

An alternative representation can be derived by observing that completely e- and 
f-bonded graphs are uniquely related 17 (e) ~ lr(/) by the substitution e = f+ 1 and 
ignoring all diagrams of lower order |f7^^| < |f7|; 

r,(e) = r,(/)+ £ r,Uf)- (7) 

k 

This relation allows to replace 17 (e) in (|^ by 17(/) and to change the representation of 
Dx from the fully connected e-bonded graphs to Mayer diagrams. The resulting opera¬ 
tor then applies to the virial expansion of Q, where the functional derivative substitutes 
any subgraph 17(/) C r„ x{f) in a Mayer diagram by 17(e) 

which reproduces the definition of the normalized r-particle functional. At the level of 
individual diagrams, this operation can also be written as 

rW(/,e) = r,(.)^^p,(/). (9) 

Examples for the 2-rooted diagrams of fourth virial order are shown in Fig. [^). 

In the following, we will use the close relationship between the virial expansion (j^ 
and the functional derivative (|^ to transfer the methods previously developed for the 
free energy to derive the r-particle correlation functionals. This approach is divided into 
two steps: the first one translates the Mayer into intersection diagrams, while the second 
determines their intersection probabilities. Let us first summarize the central ideas and 
notations from ifT^ : 

The approximation method for the free-energy functional uses an expansion in the 
number of intersection centers, whereas the virial series is an expansion in increasing 
powers of the particle density. Both are uniquely related by Ree-Hoover (RH) diagrams 
r„ X, which derive from Mayer graphs i7jt by inserting 1 = etj — fij for each pair ij 
of nodes not bonded by an f-function O18ill9ll20l . Their respective cluster integrals are 
related by the “star-content” x as introduced in ifTSl : 

k' k 


(10) 
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which satisfies the recursion relations 

7t-\r„) = {n-2)r„^i, = (11) 

under removal of a node and its associated f-bonds from a labeled diagram ^ > 

This map can be reversed for a nontrivial RH-graph by adding a node and bond¬ 
ing it by f-functions to all previous vertices % : f),_i For any lowest element 

= 0, this defines a unique RH-class 


^no.k=\J^"i^no,k), (12) 

m=0 

whose intersection networks can be contracted into a common pattern with an equal 
number of intersection centers. Each intersection diagram therefore belongs to a unique 
RH-class, whose elements can be summed up into a generic functional, weighted by the 
numerical prefactors of the virial expansion 

~ 1 _ l^«l 

— ,^n,kk^n,k-f ^n,k — i. \| j 1^-^/ 

n! |Aut(r„_<.)| 

where the symmetry factor (7„ ^ counts the number of inequivalent labelings, determined 
by the coset of permutations and the automorphism group Aut(f7; ,t). 

The approximate free-energy functional derives from the sum over all contracted 
intersection diagrams of a RH-class, where the particle domains ,... C 
of index I — intersect in the center G D,j n ...n. Its functional then 

factorizes into a convolute of integral kernels K, satisfying the rules 

K{y!,^+yi^)=K{yi^)+K{yl^). K{yi^yi^)=K{yi^)K{yi^). 

K{ya[e...e\)=K{yi)[e...e\ 

for the intersection networks yi', yj^, and the product of Boltzmann functions [e ...e\. 

The integral kernel for A-dimensional particles in the A-dimensional, Euclidean 
space is a combination of A -I- 1 weight functions determined by the func¬ 

tional derivative with respect to the weight function of the particle domain Wq 

^{Ya "') = Wq {rail ) ' ' ' ^0 i^ai„) , (15) 

where the derivative at intersection point rai = ra — rt is defined by 


N 


- L L ^! ■ 


H-ik( 


k=\ 


' ai\i’ 


■^'aik 


5wo‘(ra;,)...5wo*(r„4) 


( 16 ) 


These definitions provide a set of rules, which are sufficient to derive approximations 
of the free-energy functional for any number of intersection centers. And, as will be 
shown in the following, they also apply to the correlation functionals, replacing star- 
diagrams by rooted graphs and a corresponding change in their automorphism groups. 

The main step in the construction of the correlation functional is again the trans¬ 
formation of the virial series 0 from rooted Mayer diagrams to a corresponding set 
of intersection networks. And in analogy to the star-graphs, this transformation re¬ 
quires the intermediate step of inserting 1 = eij — ftj for each pair i, j of nodes not 
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bonded by f-functions, resulting in a change of the virial series from Mayer to r-rooted 
RH-diagrams 11911211 . Examples for 2-rooted graphs of the fourth virial order are 
shown in Fig.[^). 

Following the conventions of ca, we define the notation for rooted diagrams; 

(r) 

Definition 1. Let denote a labeled r-rooted Mayer diagram with r white and n — r 
black nodes. The rooted points are mutually e-bonded, while the black points are 1-path 
connected in the subset off-bonds. 

A black node can be removed by deleting its vertex and all associatedf-bonds : 

ri’ leaving a residual diagram, which is either a new r-rooted 

(r) (r) 

Mayer graph F _ j or a sum of disjunct diagrams with articulation points Ff ^ 


~(r) 

Definition 2. Let Fff denote a labeled r-rooted RH-diagram with r white and n — r 
black nodes. The rooted points are mutually e-bonded, while the black points are 1- 
path connected in the subset off-bonds. 

A black node without e-bonds can be removed by deleting its vertex and all associ- 

1 — ^(r) 

ated f-bonds ; Fj"/, —>■ {f^_i 0}> leaving either a new or the trivial RH-graph. 

To rewrite the virial integrals (j^ in rooted RH-diagrams, observe that the rooting- 
process only exchanges fully f-bonded subdiagrams Fr{f) by its correspond¬ 

ing e-bonded graph i7(e), whereas the Ree-Hoover transformation only operates on 
pairs of unbonded nodes. The two operations are therefore mutually exclusive and com¬ 
mute. From this follows that the same functional derivative (|^, which transforms star- 
diagrams into rooted Mayer graphs, also applies to RH-diagrams 


rr(e) 


_^_ p , = 

5Fr{f) 


(17) 


The previously derived results for unrooted RH-graphs therefore remain valid for its 
rooted forms. 

This connection can be immediately applied to rewrite the transformation between 
Mayer and RH-diagrams 


^n.k '=\.k> 

Fin _ V rin 

^n.k L-i ^n,k' ’ ^n,k 

k' 


F.k ^\.ki 


E (- 1 )' 


n.k 


l-ir 


n,k''' 


^ n^k' 


(18) 


and to express the sum over rooted RH-diagrams 


rj'-’ = y F^l] = 

n.k' 


“n,k ^n,k 


(19) 


in terms of the “root-content” 
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which satisfies an analogous recursion relation as the star-content GD under removal 
of a black node. This is readily seen by commuting with the functional derivative 
(j^ and using ( [TT] l 

=(n-2)r„W , (21) 

with the corresponding relation for RH-diagrams, where the removal of a black node 
and its « — 1 f-bonds 


= (-1)"-' E«S,f w „ 


('-) r{r) 


= («-2)r„w 

k' 


^(0 


( 22 ) 


yields the recursion relation 

= (23) 

The successive application of * traces each diagram to a unique lowest element 

1 fr') 

71^iJ = 0, which is the first element of the RH-class 


m=0 


(24) 


^(0 




defined by the inverse map n : ^ ^ F^j^, which attaches a node to the previous n — 1 

vertices by f-bonds. The set of rooted RH-diagrams therefore separates into RH-classes, 


whose root-contents can be recursively calculated by (23 i, proving 


~{r) ~(r) ~(r) 

Lemma 1. The root-content ofF , G A , with lowest element F^ , is determined by 


(r) 

^ni 


= (-l)(2)-("2'') 


(»- 2 )! Jr) 
(no-2)! "oF ■ 


(25) 


The last step in rewriting the correlation function (|^ in RH-graphs is the transition 
from labeled to unlabeled diagrams. As for the star-diagrams ( [T3| , the symmetry factor 

(r) ' -' 

(7^ ^ counts the number of inequivalent permutations of particle indices, determined by 
the coset of the permutation and automorphism group of rooted diagrams. 

~(r) '^(r) '^(r) 

Lemma 2. Let F^j^ G denote an element of the RH-class with lowest element Ff^ 
Its inequivalent labelings of the r white and n — r black nodes are permuted by the coset 
group 

Sr X Sn-r/Aut(r}^^) 

whose automorphism group factorizes into the direct product 


(26) 
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This is shown as follows: The white and black nodes are labeled independently, result¬ 
ing in the decoupling of the permutation group Sr x S„^r, proving ( [26] l. Whereas the 
automorphism group factorizes, because each of the no nodes is linked to at least one e- 
bond, while the residual n — no nodes are completely f-bonded. Any exchange of labels 
between these two groups therefore results in an inequivalent permutation, leaving the 
n — no vertices as an invariant set under relabeling, proving ©■ 

To rewrite the virial expansion (j^ in RH-diagrams, let us define the symmetry factor 
of RH-integrals 





'n,k 


(28) 


combining the root-content < [25] t and the number of inequivalent labelings 



\Sr X /-| 

|Aut(®)| 


r\{n — r)\ 1 

(«-«o)! |Aut(f„2)| ’ 


(29) 


which yields the numerical prefactor of the rooted RH-class: 

~(r) 

Corollary 1. The symmetry factor of the r-rooted RH-diagram Fff G lowest 

element r >2 is determined by 



This result follows by inserting ( |25] l, ( p9] l into (j2^. 

Up to now, no approximation has been made on the virial expansion (|^. The next 
step is therefore to simplify the integrals by restricting the number of intersection cen¬ 
ters in which the particles are allowed to overlap. A useful observation is the following 
property of intersection diagrams of a RH-class: 

^(r) 

Lemma 3. The intersection network of the class is defined by its lowest subgraph. 

The proof begins with the intersection diagram of the lowest element. Any further par¬ 
ticle, added by tt, can then be chosen to overlap with the previous intersection centers. 
This shows that a new particle can be added without changing their number. Repeated 
operation with n then completes the proof. 

The intersection diagram with the lowest number of intersection centers therefore 
defines the “backbone” diagram for the entire RH-class. This is the basic idea for the 
resummation of RH-diagrams of a given RH-class and the approximation transferring 
the virial series ^ to the generic r-particle correlation functional: 

Theorem 1. Let Ff^ ^ denote the lowest element of the RH-class The generic r- 

particle correlation functional of the intersection network ... 7 /^ [e... e] with parti¬ 
cle indices I G (ii,..., is determined by 


giy-.-ir IP [rty ,..., r,Jr„,,..., ) = F^^^ {e)5n^,j5Ai + (-1) 


I 


r\ a 


r) 

nQ,k 





( 31 ) 
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with the notation 5 ai = )... 5(ra^/^)5(ra^^j ).. .5 (rap) for the product of delta- 

functions. 


This result follows from inserting (30 1 into the cluster expansion ^ and the cancelation 
of signs due to the identity 






(32) 


which leaves an overall constant, depending only on the lowest element of the RH-class 
and an n-dependent binomial coefficient. An exception provides the leading diagram 
no = r of each r-correlation functional. Without an f-bond, its intersection probability 
vanishes and therefore has to be included separately. Finally, the delta-functions 5 ai 
have been added in ([3T| to achive a symmetric formulation of the integrals 


{8iii2fii,i2]irii^ri2) = j giui2Wi\^ri2\ra,rb,rc)fip2irh,rif)d'Yad'Yhd'Yc 

[giihfhhMra.rb^rc) = / gi^^if3{ri^,rifra,rb,rc)fip^{ri^,ri^)dyi^d'Yi^ 


(33) 

(34) 


in the particle and intersection coordinates, which proves useful in later applications. 


3. Examples of R-Particle Correlation Functionals 


This last section presents four examples. We begin with the explicit derivation of the 
2-particle correlation functional with two intersection centers and compare its contact 
probability with the Ornstein-Zernike solution of Wertheim, Baxter, and Thiele. The 2- 
and 3-particle correlations are then calculated for up to four intersection centers and 
compared to the Kirkwood superposition approximation. 

The first f-bonded RH-diagram of the 2-particle correlation functional ^ i > 

whose two white and one black nodes define a backbone diagram with two intersection 
centers. All further networks of the same RH-class can then be contracted to the pattern 



^i\i2 “f 


Yjii-inyi2h-in 


(35) 


Applying the rules (14 1 and •ED for the representation of the intersection kernel and 
introducing the notation wj, = WQ{rai) for the volume weight, the virial series sums up 
to the generating function 


giii2\2{rii,ri2\ra,rb) 

; J 


'2l3-ln 


P 13 ■ • ■Pindjii ■ ■ ■dYi„ 


n> 2 i 3 ...i„ “ 


= ei 


102 


i^AI + E E ■■■^b P'3 ■ ■ ■ Pin d7i3---d7i„) 


(36) 


n>2i2.-.in 

= ei02 i^AI + %% K' ^b E ’ 

n>2 


where we introduced the jc-variable 




...ap = J Wo(r„,,). ..Wo{rapi)pidyi 


(37) 
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from ca and used the property of ( [T6] l that the intersection probability for a single 
particle is zero ^aw'a = 0. The final correlation functional has then the analytic form 

I’l i2 

/ W/7 VV, \ 

gili2\2(ri^,ri^\ra,r|,)=ei^iA5{rai^)5{rhi^) + ^a% 2 — -) ■ (38) 

This result is exact up to the third virial order, but significantly improved by the 
additional pole at x^b = 1. To illustrate this effect of the resummation process, let us 
derive the leading order of the contact probability (|r,, — rjj | = D+) for spheres of 
radius R and diameter D — 2R. The polynomial expansion of the numerator of ( [38| up 
to first order in p yields three terms 


'1 *2 

f / tVa Wr \ 

gilh 12, ri2) = e;, i2 [Sai + -) djadjb 

J \ i Xnh' 


-lll2 


J 


e 


Sai + ^c, 
hh 


w’J J w'i &b{w'^w‘l )p,3 dji 


1 Xab' 
b "’b )P'3 ^><3 


dYadYb 


(39) 


(1 -Xab)^ 


(1 -X„fo)2 

(^5Aiil-2xab) + J{%w‘^w‘^){%w‘^w‘^)pi^dYi 2 + dYadYb , 


of which two correspond to the second and third virial orders, whereas —2xab origi¬ 
nates from resummation. Their respective integrals determine the lens-like volumes of 
overlapping spheres 


y«(r) = "^R^ 


1- 


3 r 1 


■ (- 
4R 16\R 


3n 


Voir) = yZ)' 


1 -- 


3 r 1 


■ (- 
AD 16\D 


(40) 


and can be readily evaluated for |r,j — | = D. Observing that the last term corresponds 

to the third virial integral 


J {^aW‘Jw']){^bW^W^)Pi^dYi^ =5 ai 


fiihfhhPhdYh = SaiP VD(r,j,2 = D) (41) 


and that the intersection coordinates can be chosen to coincide 


J 2Cab^ {^aii ')^iXbi2 ^dYadYb — P^Ri^ab — 0) — T] , 


(42) 


the integral (391 yields the contact probability as a function of the packing fraction 77 ; 

••2 


ghhinih D+) ^’^■*'2^^'^^((l-77)0 


(1-77)2 

1 + 2^7 

(1-7?)2 


H” ( 


( ) 


(43) 


3 ) ’ 


which agrees with the Carnahan-Starling polynomial to first order in the numerator as 
well as the Wertheim, Thiele, Baxter solution of the Ornstein-Zernike equation Kni 
[33][35][36l . The latter observation is especially interesting because its closing condition 
C 2 {r > D) =0 corresponds to the 1-center approximation of the functional expansion, 
whereas ( [4T| reflects the 2-center representation of g 2 - To obtain the same accuracy for 
the distribution functional therefore requires a larger number of intersections centers 
than for its dual direct correlation. 
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Deriving higher order functionals is now a matter of simple algebra. Here we list 
the leading orders of the 2-particle correlations with up to four intersection centers. The 
corresponding intersection diagrams for the RH-classes in the notation of Fig. [^) are 
obtained by successive contraction of pair-wise intersection centers 


A • p- ■ -X-p- ■ "Jlh-’-in nj2h---hi 

^‘2,1 • la lb 

+ (44) 

IP) , r{2) . yi'3V3'2 . Jlhh-anyiiUh-inyhhh-in hkh-in 

2^2.1+2143 . e,,,2+^!l(2la lb +^u'2/ai 102 Ibi lb 2 

_Lp. . p. . Vl'4*5-in V2'3'5-'n j2'4*5-'« 

A^erii2C(34 ffll faj lb[ lb2 ’ 


with the corresponding 2-particle correlation functionals 


gki2\2=eiJ5Ai + %%^^) 

' ^ ^ab' 

giihb = e/,12 (^21/ + ^a%w 2 ^b 2 <^ab + %^b^cWa'w‘^w‘^' ^a^ 


^abc 


i\ h iih 

(1 -XabcY-' 


gk 12 14 = e /,,'2 (5a/ + w^\ Xa,bi 


+ 9a, 9a2 %, 9b2 (< ^*2 ) 


■ ^a[a2b\b2 

+ ^9 9 9.9. <<<yb2yaM2 \ 

+ ^9a,9a2^b,Jb2 

where we introduced the Boltzmann weighted densities 


(45) 


yabc = 11 f ei,uei^i2w‘]w^ Pi^Yi^dYi^ 

ya,a2b,b2 = 11 e;3;4W^3 w;^ . 

hM'' 


(46) 


As a last example, we present the leading term of the 3-particle correlation functional 
with 3 intersection centers. Its intersection diagram has the form 


A 


• p- ■ p- . p- ■ -1- P ■ ■ P‘ ■ P‘ ■ 3X16 ■■■/«'y^ 24 -"lll 32^3/4-■■/ll 

3,1 ' ^ fa th 5 


with the corresponding correlation functional 

8i\i2h I 3 “ ^hh^hh^hh 


h o ^3 
(1 XabcY ^ 


(47) 


(48) 


The contact probability is therefore to leading order gs (1 — p) which shows that 
the Kirkwood approximation « ^/i/ 2 ^/ 2 / 3 ^/ 3 /i = [g2? -- (1 -is not appli¬ 
cable when all three particles are close together, in accordance with results obtained 
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from computer simulations ||32|. It is not difficult to generalize this result to an ar¬ 
bitrary r-particle correlation function, whose leading term is ~ (1 ~ whereas 
the Kirkwood approximation suggests ^ (1 — This excludes the 

superposition approximation as a construction principle for any distribution functional. 


4. Discussion and Conclusion 

The current derivation of the distribution functionals completes our investigation of the 
hard-particle correlations, which started with the analysis of Rosenfeld’s fundamental 
measure theory. Together with this previous result, it is now possible to calculate all 
direct and distribution functionals in an expansion of intersection kernels, highlight¬ 
ing the hard-particle interaction as the only known potential whose correlations can be 
derived by analytic methods instead of molecular dynamic or Monte Carlo simulations. 

Despite their complex structure it is possible to derive several characteristic prop¬ 
erties from these functionals. For example, we have shown that the Wertheim, Thiele, 
Baxter solution only includes the third virial integral, whereas the 2-center approxima¬ 
tion also contains contributions from the 4-particle diagram. The expansion of the pair 
correlation to leading order in the intersection kernels therefore applies not only to more 
general particle geometries than Wertheim’s solution, but also includes approximations 
of higher virial terms. We also tested the Kirkwood superposition approximation, prov¬ 
ing that it fails for any distribution functional, a result that has been previously shown 
numerically only for the 3-particle correlations. 

These examples demonstrate the difficulties to find appropriate methods to eval¬ 
uate the non-local correlation functionals. A possible approach is, e.g., to generalize 
Wertheim’s ansatz of expanding the intersection kernels in convolutes of r-point den¬ 
sities. But the integrals soon become intractable to solve. The only alternative, which 
also applies to concave geometries, is the triangulation of the surfaces and to calculate 
the functionals numerically. The intersection configurations of the particles can then be 
generated either by discretization of the imbedding space or by Monte Carlo sampling 
as has been done by Ree and Hoover. But the time-consuming step remains the determi¬ 
nation of the intersection domains and their Euler densities. Both problems have been 
solved in the framework of discrete geometry and implemented in software libraries 
for 3-dimensional image manipulation. But the vast number of intersection configu¬ 
rations that have to be evaluated while minimizing the grand potential requires further 
approximations to reduce the calculational costs in order to be competitive to molecular 
dynamic and Monte Carlo methods. 

The only particle geometry for which the functionals can be evaluated by algebraic 
relations are spheres. Instead of looking for approximations of the functionals, it is 
therefore more appropriate to replace the geometry of molecules by sets of overlap¬ 
ping balls. This approach will be discussed in a forthcoming article, where the grand 
potential for realistic molecules of hard and soft interactions will be derived ifTSll . 


Acknowledgements. The author wishes to thank Andre Bardow and Kai Leonhard for their support of this 
work. 


References 

1. Baxter, R.: Ornstein-zernike relation for a disordered fluid. Aust. J. Phys. 21, 563-569 (1968) 

2. Blaschke, W.: Vorlesungen iiber Integralgeometrie. Deutscher Verlag der Wissenschaften (1955) 


Distribution Functionals for Hard Particles in N Dimensions 


13 


3. Brader, J.M., Esztermann, A., Schmidt, M.: Colloidal rod-sphere mixtures: Fluid-fluid interfaces and the 
onsager limit. Phys. Rev. E 66(3), 031,401 (2002) 

4. Chem, S.S.: On the kinematic formula in the euclidean space of n dimensions. Am. J. Math. 74, 227-236 
(1952) 

5. Chem, S.S.: Integral formulas for hypersurfaces in euclidean space and their applications to uniqueness 
theorems. Indiana Univ. Math. J. 8, 947-955 (1959) 

6. Chem, S.S.: On the kinematic formula in integral geometry. J. Math. Mech. 16, 101-118 (1966) 

7. Evans, R.: The nature of the liquid-vapour interface and other topics in the statistical mechanics of non- 
uniform, classical fluids. Adv. Phys. 28, 143-200 (1979) 

8. Evans, R.: Fundamentals of Inhomogeneous Fluids, chap. 3, pp. 85-175. Marcel Dekker, New York 
(1992) 

9. Hansen-Goos, H., Mecke, K.: Fundamental measure theory for inhomogeneous fluids of nonspherical 
hard particles. Phys. Rev. Lett. 102(1), 018,302 (2009) 

10. Hansen-Goos, H., Mecke, K.: Tensorial density functional theory for non-spherical hard-body fluids. J. 
Phys.: Condens. Matter. 22, 364,107 (2010) 

11. Korden, S.: Deriving the rosenfeld functional from the virial expansion. Phys. Rev. E 85(4), 041,150 
( 2012 ) 

12. Korden, S.: Density functional theory for hard particles in n dimensions (2014) 

13. Korden, S.: Lattice-fluid models derived from density functional theory (2015) 

14. Marechal, M., Goetze, H.H., Hartel, A., Lowen, H.: Inhomogeneous fluids of colloidal hard dumbbells: 
Fundamental measure theory and monte carlo. J. Chem. Phys 135, 234,510 (2011) 

15. Marechal, M., Lowen, H.: Density functional theory for hard polyhedra. Phys. Rev. Lett. 110, 137,801 
(2013) 

16. Marechal, M., Zimmermann, U., Lowen, H.: Freezing of parallel hard cubes with rounded edges. J. 
Chem. Phys. 136(14), 144,506 (2012) 

17. McDonald, I.R., Hansen, J.P.: Theory of Simple Liquids. University of Cambridge (2008) 

18. Ree, F.H., Hoover, W.G.: Fifth and sixth virial coefficients for hard spheres and hard disks. J. Chem. 
Phys. 40, 939 (1964) 

19. Ree, F.H., Hoover, W.G.: Reformulation of the virial series of classical fluids. J. Chem. Phys 41, 1635 
(1964) 

20. Ree, F.H., Hoover, W.G.: Seventh virial coefficients for hard spheres and hard disks. J. Chem. Phys 46, 
4181 (1967) 

21. Ree, F.H., Keeler, N., McCarthy, S.L.: Radial distribution function of hard spheres. J. Chem. Phys 44, 
3407-3425 (1966) 

22. Riddell, R.J., Uhlenbeck, G.E.: On the theory of the virial development of the equation of the state of 
monoatomic gases. J. Chem. Phys. 21, 2056 (1953) 

23. Rosenfeld, Y.: Scaled field particle theory of the structure and the thermodynamics of isotropic hard 
particle fluids. J. Chem. Phys. 89(7), 4272^287 (1988) 

24. Rosenfeld, Y: Free-energy model for the inhomogeneous hard-sphere fluid in d dimensions: Structure 
factors for the hard-disk (d=2) mixtures in simple explicit form. Phys. Rev. A 42(10), 5978-5989 (1990) 

25. Rosenfeld, Y: Density functional theory of molecular fluids: Free-energy model for the inhomogeneous 
hard-body fluid. Phys. Rev. E 50(5), R3318-R3321 (1994) 

26. Rosenfeld, Y: Free energy model for the inhomogeneou hard-body fluid: Application of the gauss-bonnet 
theorem. Mol. Phys. 86, 637-647 (1995) 

27. Rosenfeld, Y, Schmidt, M., Lowen, H., Tarazona, P: Dimensional crossover and the freezing transition 
in density functional theory. J. Phys.: Condens. Matter 8, L577-L581 (1996) 

28. Rosenfeld, Y, Schmidt, M., Lowen, H., Tarazona, R: Fundamental-measure free-energy density func¬ 
tional for hard-spheres: Dimensional crossover and freezing. Phys. Rev. E 55, 4245^263 (1997) 

29. Santalo, L.A.: Integral Geometry and Geometric Probability. Addison-Wesley (1976) 

30. Schmidt, M.: Fluid structure from density-functional theory. Phys. Rev. E 62(4), 4976^981 (2000) 

31. Tarazona, P., Rosenfeld, Y: From zero-dimensional cavities to free-energy functionals for hard disks and 
hard spheres. Phys. Rev. E 55, R4873-R4876 (1997) 

32. Taylor, M.P., Lipson, J.: On the born-green-yvon equation and triplet distributions for hard spheres. J. 
Chem. Phys. 97, 4301^308 (1992) 

33. Thiele, E.: Equation of state for hard spheres. J. Chem. Phys. 39, 474-479 (1963) 

34. Uhlenbeck, G.E., Ford, G.W.: The Theory of Linear Graphs with Applications to the Theory of the Virial 
Development of the Properties of Gases, Studies in Statistical Mechanics, vol. 1. Interscience Publishers 
(1962) 

35. Wertheim, M.S.: Exact solution of the percus-yevick integral equation for hard spheres. Phys. Rev. Lett. 
10, 321-323 (1963) 

36. Wertheim, M.S.: Analytic solution of the percus-yevick equation. J. Math. Phys. 5, 643-651 (1964) 

37. Wertheim, M.S.: Fluids of hard convex molecules i. basic theory. Mol. Phys. 83, 519-537 (1994) 



14 


Stephan Korden 


38. Wertheim, M.S.: Fluids of hard convex molecules ii. two-point measures. Mol. Phys. 89, 989-1004 
(1996) 

39. Wertheim, M.S.: Fluids of hard convex molecules iii. the third virial coefficient. Mol. Phys. 89, 1005- 
1017(1996) 

40. Wertheim, M.S.: Third virial coefficient of hard spheroids. Mol. Phys. 99, 187-196 (2001) 



